An integrative model of motor unit activity during sustained submaximal contractions.

The purpose of the study was to expand a model of motor unit recruitment and rate coding (30) to simulate the adjustments that occur during a fatiguing contraction. The major new components of the model were the introduction of time-varying parameters for motor unit twitch force, recruitment, discharge rate, and discharge variability, and a control algorithm that estimates the net excitation needed by the motoneuron pool to maintain a prescribed target force. The fatigue-induced changes in motor unit activity in the expanded model are a function of changes in the metabolite concentrations that were computed with a compartment model of the intra- and extracellular spaces. The model was validated by comparing the simulation results with data available from the literature and experimentally recorded in the present study during isometric contractions of the first dorsal interosseus muscle. The output of the model was able to replicate a number of experimental findings, including the time to task failure for a range of target forces, the changes in motor unit discharge rates, the skewness and kurtosis of the interspike interval distributions, discharge variability, and the discharge characteristics of newly recruited motor units. The model output provides an integrative perspective of the adjustments during fatiguing contractions that are difficult to measure experimentally.